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We calculate the velocity of sound in an ideal gas of massive hadrons with non-vanishing 
baryon number. The gas is in thermal and chemical equilibrium. Also we show that the 

temperature dependence T(r) = To ■ (2M a is approximately valid, when the gas expands 
longitudinally according to the Bjorken law. 



I. INTRODUCTION 

The main subject of the following paper is the calculation of the sound velocity c s in the multi-component 
hadron gas with nonzero baryon number density. The sound velocity appears in the hydrodynamical description 
of the matter which is created in the central rapidity region (CRR) of a heavy-ion collision. It turned out jjJD 
that the evolution of the matter (whatever it is: a quark-gluon plasma or the hadron gas) proceeds as longitudinal 
and superimposed transverse expansions. And the latter has the form of the rarefaction wave moving radially 
inward with the velocity of sound. The second place where the sound velocity appears is the temperature 
equation for an ideal baryonless gas (it can be a gas of quarks and gluons as well), which cools as the result 
only of the longitudinal expansion (see e.g. 0) 

T(t) = T .(^) C ' , (1) 

where r is a local proper time, tq an initial moment of time and c s is assumed constant. In the following, we 
show that the above equation, now as some approximation, is also valid in some range of temperature for the 
hadron gas with non- vanishing baryon number. 

In our previous paper Q we showed that taking into account also the transverse expansion (in the the form 
of the rarefaction wave) changes the J/^> theoretical suppression pattern qualitatively. This was done for the 
scenario without the quark-gluon plasma appearance in the CRR. Of course, this is very important because 
J/'f suppression observed in NA38 H and NA50 |(| heavy-ion collision data is treated as the main signature of 
the creation of the quark-gluon plasma (QGP) during these collisions 0. Since the study of J/^> suppression 
patterns was the main purpose of our paper [0 , we left aside the more complete presentation of the behaviour 
of the sound velocity in the hadron gas. Now, we would like to present our results in the form of a separate 
paper. 



II. THE EXPANDING MULTI-COMPONENT HADRON GAS 

For an ideal hadron gas in thermal and chemical equilibrium, which consists of I species of particles, energy 
density e, baryon number density Ub, strangeness density ns, entropy density s and pressure P read (h = c = 1 
always) 



2 - 2 hC Jo exp{^} 
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where Ei = (m 2 + p 2 ) 1 / 2 and m„ B i: Si, p~i, Si and gi are the mass, baryon number, strangeness, chemical 
potential, spin and a statistical factor of specie i respectively (we treat an antiparticle as a different specie). 
And pn = Bi^iB + Si /is j where fiB and us are overall baryon number and strangeness chemical potentials 
respectively. 

We shall work here within the usual timetable of the events in the CRR of a given ion collision (for more 
details see e.g. II). We fix t = at the moment of the maximal overlap of the nuclei and assume that it takes 
place at z = 0, where z is the coordinate of a collision axis. After half of the time the nuclei need to cross each 
other, matter appears in the CRR. We assume that soon thereafter matter thermalizes and this moment, tq 
(t = y/t 2 — z 2 ) , is estimated at about 1 fm . Then matter starts to expand and cool and after reaching the 
freeze-out temperature I/. . it is no longer a thermodynamical system. As we have already mentioned in the 
introduction, we assume that this matter is the hadron gas, which consists of all hadrons up to Q~ baryon. The 
expansion proceeds according to the relativistic hydrodynamics equations and for the longitudinal component 
we have the following solution (for details see e.g. 

/ x s o T o i \ n° B T 

s(t) = , n B {T) = -2— , (3) 

r r 

where sq and n° B are initial densities of the entropy and the baryon number respectively. The transverse 
expansion has the form of the rarefaction wave moving radially inward with a sound velocity c s 0,0 . 

To obtain the time dependence of temperature and baryon number and strangeness chemical potentials one 



has to solve numerically equations (2b - 2d) with s, ns and n$ given as time dependent quantities. For s(r), 
ub{t) we have expressions (^) and n$ = since we put the overall strangeness equal to zero during all the 
evolution (for more details see |^]). 

III. THE SOUND VELOCITY IN THE MULTI-COMPONENT HADRON GAS 

In the hadron gas the sound velocity squared is given by the standard expression 

Since the experimental data for heavy-ion collisions suggests that the baryon number density is non-zero in the 
CRR at AGS and SPS energies [^0|-^2|, we calculate the above derivative for various values of iib- 

To estimate initial baryon number density n B we can use experimental results for S-S jl0| or Au-Au |Il| , |l2"| 
collisions. In the first approximation we can assume that the baryon multiplicity per unit rapidity in the CRR is 



proportional to the number of participating nucleons. For a sulphur-sulphur collision we have dNs/dy ^6 |10 



and 64 participating nucleons. For the central collision of lead nuclei we can estimate the number of participating 
nucleons at 2A — 416, so we have dNs/dy = 39. Having taken the initial volume in the CRR equal to ~kB? a ■ 1 fm, 
we arrive at n B ~ 0.25 /to -3 . This is some underestimation because the S-S collisions were at a beam energy 
of 200 GeV/nucleon, but Pb-Pb at 158 GeVynucleon. From the Au-Au data extrapolation one can estimate 
n B = 0.65 jmT 3 These values are for central collisions. So, we estimate (|]) for tib = 0.25, 0.65 fm~ 3 and 
additionally, to investigate the dependence on tib much carefully, for ns = 0.05 fm~ 3 . The results of numerical 
evaluations of (Q) are presented in Fig. |l|. For comparisomwe drew also two additional curves: for jib = and 
for a pure massive pion gas. These curves are taken from 




FIG. 1. Dependence of the sound velocity squared on temperature for various values of ns: ns = 0.65 fm 
(short-dashed), tig = 0.25 fm~ 3 (dashed), ns = 0.05 fm~ 3 (solid) and ns = (long-dashed). The case of the 
pure pion gas (long-long-dashed) is also presented. 



The "physical region" lies between more or less 100 and 200 MeV on the temperature axis. This is because 
the critical temperature for the possible QGP-hadronic matter transition is of the order of 200 MeV jl3| and 
the freeze-out temperature should not be lower than 100 MeV pT] , In low temperatures we can see completely 
different behaviours of cases with n B = and n B ^ 0. We think that this is caused by the fact that for n B ^ 
the gas density can not reach zero when T — » 0, whereas for n B — it can. For the higher temperatures all 
curves excluding the pion case behave in the same way qualitatively. From T w 70 M eV they decrease to their 
minima (for n B = 0.65 /m" 3 at T S 221.3 MeV, for n B = 0.25 /m" 3 at T S 201.5 MeV, for n B = 0.05 /m" 3 
at T = 176.6 MeV and for n B — at T = 160 MeV) and then they increase to cover each other above 
T ps 250 MeV. To investigate this problem we check the content of the hadron gas when the temperature 
changes. 
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FIG. 2. Fraction of mesons in the multi-component 
hadron g function of temperature for 

ub = 0.65 fm' 3 (short-dashed), ns = 0.25 fm~ 3 
(dashed), tlb = 0.05 fm' 3 (solid) and n_g = 
(long-dashed). 
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FIG. 3. Fraction of pions in the multi-component 
hadron g ELS cLS cL function of temperature for 
ub = 0.65 fm~ 3 (short-dashed), ns = 0.25 fm~ 3 
(dashed), ns = 0.05 fm~ 3 (solid) and ns = 
(long-dashed). 
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FIG. 4. Fraction of rj, p and U) mesons in the 
multi-component hadron gas as a function of tem- 
perature for n B = 0.65 /m~ 3 (short-dashed), 
ns = 0.25 /m~ 3 (dashed), n B = 0.05 fm' 3 (solid) 
and n B = (long-dashed). 
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FIG. 5. Fraction of strange particles in the 
multi-component hadron g £LS EIS ct function of tem- 
perature for n B = 0.65 fm~ 3 (short-dashed), 
n B = 0.25 fm~ 3 (dashed), n B = 0.05 fm~ 3 (solid). 
The curve for n B = (long-dashed) covers the curve 
for n B — 0.05 fm~ 3 almost completely. 



We define the fraction of some kind of particles in the hadron gas as 

U = - . (5) 
Q 

where g p is the density of these particles and g the density of the hadron gas. We calculate the above-mentioned 
fractions of: a) mesons; b) pions; c) 77, p and uj mesons and d) strange particles. The results are presented in 
Figs. §§. 

Coming back to the behaviour of the sound velocity of the hadron gas above T « 70 MeV (see Fig. [lj), we 
can state that the appearance of heavier particles, especially 77, p and lu mesons, is responsible for it. First of 
all, the curve for the pure pion gas in Fig. |l| is increasing everywhere, so we have to take into account also 
other species to obtain the sound velocity decreasing in some range of temperature above T ss 70 MeV . If we 
assume that baryons cause the mentioned behaviour of the sound velocity, it would not agree with the case of 
ns = 0, because the baryon fraction is increasing with the temperature everywhere for this case, see Fig. || 
up-side-down. The same would happen for strange particles (with the only difference that the strange particle 
fraction is increasing with the temperature for all cases of ns, see Fig. ||). Only the fraction of ?y, p and ui 
mesons changes its monotonicity above T ss 70 MeV and at the temperatures similar to those where the sound 
velocity curves have their local minima (cf. Fig. and Fig. |), for all cases of n_e- 



IV. THE PATTERN OF COOLING AND ITS CONNECTION WITH THE SOUND VELOCITY 

In Sect. |n| we have explained how to obtain the time dependence of the temperature of the longitudinally 
expanding hadron gas. This dependence proved to be very well approximated by the expression 

T(r) S To ■ r~ a . (6) 

The above approximation is valid in the temperature ranges [TV ,,Tb], where T/. Q . > 100 MeV, Tq < To^ max 
and T , max » 225 MeV. We started from T equal to 221.8 MeV (for n° B = 0.65 /m~ 3 ), 226 MeV ' (for 
n% = 6.25 fm- 3 ) and 226.7 MeV (for n° B = 0.05 /m" 3 ). These values correspond to e = 3.7 GeV/fm 3 , 
which is the initial energy density in the CRR slightly above the value eo — 3.5 GeV/ fm 3 estimated by NA50 
||. Then we took several decreasing values of Tq < Tb, max - For every Tq chosen we repeat the procedure of 
obtaining the approximation (0), i.e. the power a. The values of a as a function of To are depicted in Fig. || 
for the case of ub = 0.25 fm , together with the corresponding sound velocity curve. 
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FIG. 6. The power a (crosses) from the approximation (g) and the sound velocity curve for ns = 0.25 fm 3 . 

We can see that in the above-mentioned interval of temperature the power a has the meaning of the sound 
velocity within quite well accuracy. We would like to add that we arrived at the same conclusion also for cases 
with n,B = 0.65 /m~ 3 and 115 = 0.05 /to -3 . Note that the higher temperature, the worse accuracy. This results 
from the fact that for a higher temperature the wider interval is used to obtain the formula (g). Therefore the 
approximation is worse. 

We can formulate the following conclusion: in the "physical region" of temperature and for realistic baryon 
number densities, the longitudinal expansion given by (|3|) results in the cooling of the hadron gas described by 

(g) with a = c 2 .(T a ). Note that Tjr) = T • (^) Cs is the exact expression for a baryonless gas with the sound 
velocity constant (for details see |2|,||). 



V. CONCLUSIONS 



In this paper we have calculated the sound velocity in the multi-component hadron gas with non-zero baryon 
number density. Such a gas, instead of the QGP, could have appeared in the CRR of heavy-ion collisions at 
AGS and SPS energies. We have compared the results with the sound velocities in a gas with ub = and in a 
pion gas. We have shown that in the "physical region" (temperatures between roughly 100 and 200 MeV) the 
sound velocities for various cases of the multi-component hadron gas behave qualitatively the same: firstly they 
decrease and then they start to increase with the temperature. Comparison with the case of the pion gas and 
analysis of fractions of particles in the multi-component hadron gas suggest that the appearance of 77, p, and ui 
mesons is responsible for that. 

The second result shown here is that for the multi-component hadron gas the cooling imposed by the longi- 
tudinal expansion can be described by the approximation 

r(r)ST .(J) (7) 

where Tq belongs to " physical region" . The approximation is valid for the gas with the non-zero baryon number 
density. The same conclusion was drawn for the case with jib — in || . 

The values of the sound velocity obtained in the "physical region" and the validity of approximation (Q) 
therein, cause that the cooling of the hadron gas is much slower than the cooling of an ideal massless gas with 
iib = (where c 2 s = -|). This fact has the straightforward consequence for J/* suppression: the longer hadrons 
last as a gas, the deeper suppression takes place (for details see 0). 
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